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Abstract 

We estimate the rates for radiative transitions of the lightest scalar mesons /o(980) 
and oo(980) to the vector mesons p and oj. We argue that measurements of the 
radiative decays of those scalar mesons can provide important new information on 
their structure. 


1 Introduction 


Although studied since many decades, the lightest scalar mesons and, es¬ 
pecially the / 0 (980) and a 0 (980), are still subject of debate regarding their 
fundamental structure. For example, these two mesons could be viewed as 
natural candidates for the genuine 1 3 P 0 QQ states predicted by the standard 
quark models [1], However, due to the proximity of these states to the KK 
threshold, a significant if not dominant qqqq configuration is expected from a 
phenomenological point of view. Thus, it was argued (see, e.g., Ref. [2]) that 
the genuine qq 3 Pq nonet could be somewhere near 1.5 GeV, while the states 
around 1 GeV are due to a strong S'-wave attraction between the two quarks 
and two antiquarks. 

In such a scenario the / 0 (980) and a 0 (980) mesons could be realized either in 
form of compact qqqq states [3,4] or in form of loosely bound KK states. To 
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complicate things further, within the latter picture there are even two possi- 
bilites regarding the nature of those scalar resonances, which are connected 
with the question whether there are sufficiently strong f-channel forces so 
that KK molecules are formed, as advocated in Refs. [5,6,7,8], or whether 
the meson-meson interaction is dominated by s-channel states. In the former 
case the /o(980) and do (980) would be purely composite particles, whereas 
in the latter case they would contain both elementary states and composite- 
particle components. For a much more thorough discussion on that issue and 
an overview of the extensive literature we refer the reader to the reviews in 
Refs. [9,10,11,12,13,14], 


Over the years many experiments have been proposed in order to distinguish 
among those scenarios, however, so far the smoking gun experiment has not 
been identified yet. For example, about a decade ago it was believed that 
data on the decays J/\ 1/ — > (pnn/^KK would allow to resolve the puzzle of 
the scalar mesons [15]. Even earlier Achasov and Ivanchenko had argued that 
measurements of the radiative decays of the 0(1020) to scalar mesons would 
provide decisive information on the structure of these long-debated objects 
[16]. The authors of that work demonstrated that the spectrum of, e.g., 7r°n° in 
the reaction 0 —> yA —» y7r 0 7r° would look drastically different in the presence 
or absence of a significant KK contribution, due to the proximity of both the 
mass of the scalar meson S and that of the 0 to the KK threshold. And indeed, 
the data [17,18,19] unambiguously show a prominent KK contribution. Based 
on large -Nc considerations, this was interpreted then as a proof for a compact 
four-quark nature of the scalar mesons [20]. However, it is not clear a priori 
how quantitative the large-A^c counting rules are in the scalar sector — for 
example, the large-A^c analysis of the unitarised chiral perturbation theory 
amplitudes leads to large uncertainties for the a 0 //o(980) states [21], After 
all, it is well known that, for the scalars, due to the presence of the nearby 
strong A-wave KK channel, unitarity corrections are large, as seen from the 
corresponding Flatte distributions [22,23]. Thus, the large- Nq picture might 
be obliterated. In particular, the admixture of a KK component in the scalar 
wave function should be large, as discussed in Refs. [5,6,24,25]. 


For obvious reasons, a dominant role of the KK component in the 0 —> yA 
decays is naturally expected in the scenario where the scalar mesons are KK 
molecules. Thus, it might be not too surprising that explicit calculations uti¬ 
lizing such molecular models [26,27,8] were able to describe the spectrum of 
the radiative 0 decays. Indeed, it was shown [28] recently on rather general 
grounds that, contrary to earlier claims [29,30], the available experimental in¬ 
formation is completely consistent with the molecular structure of the scalar 
mesons. We thus conclude that the radiative 0 decays measure the molecular 
component of the scalar mesons. However, other observables are to be found 
that allow one to understand how much compact structure there is in addition. 
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There is a drawback of the radiative 0 decays: beyond the prominence of 
the kaon loops, no further model-independent quantitative conclusion on the 
scalar mesons is possible because of the limited phase space available for 
these decays. In addition, gauge invariance forces the spectrum, for large pseu¬ 
doscalar invariant masses, to behave as cu 3 , where lv denotes the photon energy. 
With the photon energy being just around 40 MeV, only a small fraction of 
the spectral functions of the scalar mesons is visible in these reactions. As 
discussed in detail in Ref. [31], this causes uncertainties in the attempts to 
define the coupling constants and pole positions of the scalars. 

In view of the difficulties outlined above, with the present paper, we would 
like to draw attention to another class of radiative decays — namely, to the 
radiative decays of the scalar mesons themselves. In particular, we want to 
provide evidence for the following properties of the reactions S —> 7 V, where 
S denotes the scalar mesons a 0 or / 0 and V stands for the vector meson p or 
u>: 

( 1 ) both quark loops and meson loops can be of equal importance; 

( 2 ) there is significant phase space available for the final state; 

(3) since there is a sensitivity to the nonstrange contribution of the wave 
functions, a combined analysis of the 0 radiative decays, as well as those 
of the scalars should help to map out the underlying quark structure of 
the latter. 

Among those the first point is specifically interesting. It implies that, if the 
scalar mesons were predominantly qq or qqqq states, then quark loops as 
well as meson loops should yield sizable contributions to the decay ampli¬ 
tude and, consequently, they should have a significantly larger decay rate to 
vector mesons as compared to KK molecules, where only meson (KK) loops 
are present. 

In this context we also consider the decay S —> 77 and show that, in this case, 
there is again a striking difference in the reaction mechanism in the sense that 
now the quark loops dominate while the meson loops are suppressed. As a 
consequence, there is a certain pattern or hierarchy in the studied radiative 
decay reactions involving scalar mesons (0 —> 7 S, S —> 7 V, S —> 77 ) with 
characteristic differences for a compact (qq or qqqq) or molecular structure of 
those scalars. This suggests that a combined analysis of such decays, within 
a specific scenario of the scalar mesons, is actually a much more conclusive 
method to discriminate between these scenarios than just considering a single 
decay mode, like 0 —> 7 S', as it happened in the past. 

The paper is structured in the following way: In Sect. 2, we provide general 
expressions for the vertex function involving a photon, a scalar (S), and a 
vector (V) meson and for the total width of the transition S —> 7 V. In Sect. 
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3, we consider different decay mechanisms for the scalar mesons, he., qq and 
qqqq quark loops and meson loops, and we evaluate the decay width within 
the corresponding transition mechanisms. Sect. 4 is devoted to the decay of 
the scalar mesons to the 77 channel. Our results are analysed and discussed 
thoroughly in Sect. 5. The paper ends with a brief summary. 


2 Some generalities 


In order to write down the effective vertex for the SV 7 coupling, one is to 
respect gauge invariance for the photon. This is most easily implemented by 
using the held strength tensor for the latter. Therefore, the most general struc¬ 
ture of the SV 7 vertex is 1 

iW = M(p 2 ,q 2 )[(k ■ e v *)(q • e 7 *) - (e v * • e 7 *)(/c • q)], (1) 

where ejj and e 7 are the polarisation vectors of the vector meson and photon, 
qn and £7 are their four-momenta, respectively, and is the scalar four- 
momentum. For the (j) radiative decays, the decay amplitude exhibits a strong 
p 2 -dependence, due to the proximity of the KK threshold to both the q f> 
mass as well as to the nominal mass of the scalar meson. However, for these 
decays, we have q 2 = m|, and thus it is not possible to investigate the q 2 - 
dependence. On the other hand, for the decays S —> 7 H and the case of kaon 
loop contributions, M(p 2 ,q 2 ) shows a significant dependence on both p 2 and 
q 2 : due to the proximity of the KK threshold to the mass of the a 0 //o mesons 
and due to the finite width of the vector mesons, especially of the p-meson. 

For stable scalar and vector mesons one could directly deduce the expression 
for the total width of the transition S —» 7 H from Eq. ( 1 ), which would be 
given by 



with my and ms being the nominal masses of the vector and scalar mesons. 
For the calculation of observables, in addition to the matrix element M, two 
more ingredients are relevant — namely, the propagator of the scalar meson 
Ds(p 2 ) and that of the vector meson, Dy(q 2 ). The latter modifies the invariant 
mass spectrum of the final state, cf. the detailed discussion in the appendix. 

The finite width of the scalar mesons makes one study the decay rates as a 
function of the invariant mass of the decaying system. Consequently, in the 
total decay width, D s (m 2 s ) appears as a weight factor. For this distribution, 

1 Here the standard normalisation of the invariant amplitude is used, like, e.g., in 
Ref [29], 
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Fig. 1. Quark loop contribution to the radiative decay. 


one would need to use parametrizations given in the literature. Note that 
using such parametrizations one could run into complications connected to 
possible interference effects between the / 0 (980) and the broad I = 0 nn 
component usually referred to as “cr” [32]. In what follows we do not consider 
this possibility and give estimates for stable vectors and scalars. 

As mentioned before, it is the transition matrix element M which is the quan¬ 
tity of interest, and we investigate it now in more detail for various scenarios. 


3 The transition matrix element M 


In this section, we discuss the properties of M in various models for the scalar 
mesons and for various mechanisms of the radiative decay. 


3.1 Contribution of quark loops 


The simplest assignment for the ao(980)//o(980) mesons is the bound qq 3 Pq 
state [1]. Correspondingly, the radiative decay proceeds via a quark loop, as 
displayed at Fig. 1. If confinement is modeled by a quark-antiquark inter¬ 
action, then the ingredients needed to calculate the transition matrix ele¬ 
ment are: i) the meson-quark-antiquark vertices, ii) the dressed propagators 
of quarks, and iii) the dressed photon-quark-quark vertex. Only if the un¬ 
derlying quark model provides these ingredients in a selfconsistent way, then 
the electromagnetic transition vertex is compatible with e.m. gauge invariance 
and the SV 7 transition amplitude takes the form of Eq. (1). Reliable calcu¬ 
lations of the quark loop contributions can be done in the framework of the 
nonrelativistic quark model. The radiative transition 3 S\(qq) <->• 3 Po(qq) is an 
El transition, and the current in the rest frame of the initial meson i, in the 
lowest approximation, is 




+ (<? 


?)• 


( 3 ) 
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The expression for the matrix element M, extracted from Eq. (3), reads 

M = \ e (Q) r if, ( 4 ) 

where 

Q=\(Qq~Qq) ( 5 ) 

is the quark charge operator, and the radial part of the dipole matrix element 
between the initial and final states reads 

rif = j r 2 drR^(r)rRi(r), (6) 

with Rij{r ) being the radial wave functions for the initial and the final states, 
respectively. 

Generally, the decay rate for El transitions between the 3 S'i and 3 Pj states is 
given by (see, e.g., Ref. [33]) 

r = WJ+E a{Q)W > t ' (7) 

for the 3 S'i Pj 7 decays, and by 

r = (8) 

for the 3 Pj S '17 decays. Here u stands for the photon energy and the 
charge factor is readily calculated for a given flavour of the initial and final 
states (n denotes u and/or d quark): 


(Q ) 2 


P, for nn —> nh with the same isospin, 
( 4, for nn —> nn with different isospins, 
|, for ss —> ss. 


(9) 


One might question the applicability of the nonrelativistic or the naively rela¬ 
tivised quark model to the mentioned decays. Nevertheless, experimental data 
can be used to estimate the needed matrix element. We may use the known 
radiative decay rate of the bona fide quarkonium /i(1285) [34], as a genuine 
3 Pi qq state made of light quarks: 

r(/i(1285)) -> yp) = 1320 ± 312keV. (10) 

As shown in Ref. [35], nonrelativistic quark models with standard parameters 
yield results for the decay /i(1285) —■> 7 p that are in good agreement with the 
data. To relate this matrix element to the ones of interest we assume SU( 6 ) 
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symmetry for the wave functions which is expected to provide a reasonable 
order-of-magnitude estimate for the rates. In this case the values of the matrix 
elements r\f are to be equal for all members of the P-multiplet. 

Then one gets from Eqs. (T)-(10): 

T(a 0 -> yen) = T(f 0 (nn) -> yp) = 125 keV, 

T(a 0 -*• yp) = T(/ 0 (nn) -> yw) = 14 keV, 
r(/ 0 (ss) —>► yen) = 31 keV x sin 2 9, 
r(/o(ss) -*• yp) = 0 , 

where sin 9 measures the (small) <j> — u mixing. 

While the expressions (7) and ( 8 ) take apparently nonrelativistic form, rel¬ 
ativistic corrections are actually included in these dipole formulae, provided 
the masses and the wave functions of the initial and final mesonic states are 
taken to be solutions of a quark-model Hamiltonian with relativistic correc¬ 
tions taken into account (see Refs. [36,37] for a detailed discussion). With 
relativistic corrections to the wavefunctions taken into account the values of 
r t j for 3 P 0 and 3f states are not equal to each other anymore. So the estimates 
( 11 ) are to be considered as order-of-magnitude ones. 

Similarly we obtain for <f> decay: 

T(0 —> yao) = 0.37 keV x sin 2 9 , 

r(0 —> y/o( s s)) = 0.18 keV, (12) 

T(</>—► y/o(nh)) = 0.04 keV x sin 2 9. 

In this context let us mention that the pure ss assignment for f 0 seems im¬ 
plausible as it implies an OZI suppression of the titi mode in the (f) radiative 
decay, so that some mixing with an nn isoscalar state is needed to reproduce 
the branching fraction of the /o(980) to 7 T 7 t. 

We would like to point out that the predictions for the decay width ratios of 
rates, 


r[(/o(nn) -> yp 0 ) : (a 0 -> ycu) : (/ 0 (nn) -»• yc o) : (a 0 -»• yp 0 )] = 

9 : 9 : 1 : 1, 


(13) 


which are readily deduced from Eq. (11), are based only on the isospin relations 
(9) and are, therefore, robust. 
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Fig. 2. Annihilation contributions to the radiative decay of qqqq scalars. 
3.2 Contributions of annihilation graphs for qqqq scalars 


In the diquark-antidiquark model [3,4] both f 0 and a o can be identified with 
snsn states belonging to a cryptoexotic 3 <E> 3 flavour nonet, 


fo 

a o{h — 0 ) 


^([su][su] + M[sd]), 

^(HN - NN). 


(14) 


The radiative decays of the qqqq states proceed via annihilation of a qq pair, 
as shown at Fig. 2. Thus, in the transition S —» yp/iv, the ss pair annihilates, 
so that one has 


r(/o -» 7 P) 
r(/ 0 -»• 7^) 


(15) 


while the <f> — u mixing could generate small nonzero values of the ratio 


F(qq ->• 7^) 

r(a 0 -> 7 p) 


~ sin 2 9. 


(16) 


On the contrary, the decay cp —> jS with the four-quark scalars (14) proceeds 
via creation of a nn pair (if one neglects the <f> — u mixing), yielding 


7«o) = 9 
r (0 -»■ 7/0) 

Note that the experimental value for this ratio is around 1/6 [35]. 


(17) 


The assumption (14) is compatible with the a 0 //o mass degeneracy. However, 
with the snsn assignment for / 0 , a superallowed decay to 7T7 t is impossible, so 
that one is forced to assume a mixing of the isoscalar snsn state with a cr-like 
nnnn state (see [4]). Note that there is no such problem for the superallowed 
decay ao —> 7 nj, since the 7 contains an admixture of the strange quark pair. 


There are no theoretical estimates of the absolute values for the qqqq radiative 
decay rates. Moreover, since no single four-quark state is unambiguously iden¬ 
tified, there is also no experimental anchor at our disposal, similar to Eq. (10), 
which could allow one to predict absolute values of these rates. 




Fig. 3. Meson loop contributions to the radiative decay of scalar mesons. Diagram 
c) is required to provide the overall gauge invariance of the amplitude. 

3.3 Contribution of meson loops 


The contribution of meson loops is shown diagrammatically in Fig. 3, where 
the diagrams a) and b ) correspond to the coupling of the photon to the charge 
of the intermediate pseudoscalar meson and the diagram c) stems from gaug¬ 
ing the decay vertex of the vector meson to two pseudoscalars. The explicit 
expressions for the corresponding matrix elements read 


W a = W b 


r d A l e^-(p + q-2l) e v *-(2l + q) 

e9s9v J (2tt) 4 ((p - l) 2 - m?p){{q - l) 2 - m 2 P ){l 2 


W c 


-2eg s gv(C*-e v *) 


d 4 l _1 

(27r) 4 ((p — l) 2 — m 2 P )(l 2 


m 2 p ) ’ 


m 2 P )' 


Adding these three we get for the amplitude M introduced in Eq. (1): 

M(my,m 2 s ) = ^ 2^2 ( 18 ) 

m 2 m 2 

where a = b = —f-, with mp being the mass of the pseudoscalar; g$ and 

Tin p Tin p 

gv are the SP + P~ and VP + P~ coupling constants, while /(a, b) is the loop 
integral function. An analytical expression for this function can be found, e.g., 
in Refs. [16,29]. The dependence of (a — b) 2 \I(a, b)\ 2 on the mass of the scalar 
meson is shown in Fig. 4. 


What pseudoscalars can be responsible for the transitions under consideration? 
The Oo meson is known to couple to 7iT/ and KK, while p, u, and (f do not 
couple to Tig. Thus, for the oo, only the kaon loop is relevant. The /o meson 
couples to 7T7T and KK, whereas the only vector meson coupling to 7r7r is the p 
meson. Therefore, for the / 0 , the pion loop could contribute together with the 
kaon loop. Nevertheless, the loop integral depends drastically on the relation 
between the initial and final meson masses and on the pseudoscalar threshold. 
For both cases of the 0 —> n a o/fo an d a 0 //o —^ 7 p/u> decays the contribution 
from the pion loop is small (see Fig. 4). Thus, in what follows, only the kaon 
loop mechanism is considered. 





Fig. 4. Dependence of the function (a — b) 2 \I(a, h)\ 2 on the mass of the scalar meson 
(in GeV) for the kaon loop (solid line) and for the pion loop (dashed line). Here 
my = rritj) for the left plot (cj) decay), and my = m p / w for the right plot (scalar 
decay). 


The only input needed to evaluate the kaon loop contribution to the radiative 
decays are the effective coupling constants gs and gy. The decay constant for 
0 —> K + K~ is readily calculated from the 0 width, 

« 1.77, (19) 

47T 

and the decay constants for p/co —> K + K~ can be estimated from that for the 
p — > 7T7T decay with the help of SU( 3) symmetry considerations yielding 

1 a 2 

9v = g P K+K- = gu>K+K- = 2 9p™ ~ 2 - 13 , ~ 0 . 36 . ( 20 ) 


The last missing ingredient is gs- In Ref. [28] 2 the value of 

7 ^- = 16 mK\JsmK = 0.6 GeV 2 (21) 

47r 

was estimated. To come to this number the mass of 980 MeV for both scalars 
was used which corresponds to a binding energy of e = 10 MeV. The quoted 
estimate is based on assuming a stable molecule formed by a pointlike inter¬ 
action in the KK channel and thus should be viewed as qualitative. Such a 
value of gs lies within the range given by various parametrisations of the a 0 //o 
propagators existing in the literature — see Tables in Refs. [22,25]. 

Within the qq and the 4-quark pictures, relations between the effective cou¬ 
plings of the a 0 and f 0 to KK can be derived readily. In the qq picture, the 

2 The value given in Ref. [28] should be decreased by a factor of 2 since only the 
charged kaons contribute to the loop mechanism of relevance here. 
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scalar mesons are 3 Pq states and, for the flavour-independent strong interac¬ 
tion, one should have, approximately, 


9a 0 KK, for fo{nn) 

9 foKK = < (22) 

[ v2 g aoK K, for fo(ss), 

though one should keep in mind that effects like the instanton-induced forces 
or an admixture of the scalar glueball in the wave function of the /o may 
destroy these equalities. 

In the four-quark model, the scalar decays are superallowed and one has, for 
the ao and /o with the quark content given by Eq. (14), the relation 

9f 0 KI\ = 9a 0 KK, (23) 

which may be distorted by the aforementioned mixing of the isoscalar snsn 
with the (j-likc state. 


On the other hand, estimates for the absolute values of the scalar coupling 
constants gg involve calculations of strong decays of quark-antiquark, four- 
quark or molecular states, which are very model dependent. One might con¬ 
sider to rely on experimental data for determining the absolute values of gg. 
However, the coupling constants extracted from data are afflicted with large 
uncertainties, he., they exhibit variations up to a factor of 2-=- 3 — see Tables 
in Refs. [22,25]. This is due to the scaling property of the Flatte distributions 
near the KK threshold, as discussed in detail in Ref. [22], In the following we 
use the value of Eq. (21) for the scalar coupling. 

With the given values for the couplings one obtains, in the kaon loop model, 
the values 


T( 0 - 

■» 7 <S) = 0.6 keV, 

(24) 

r(a 0 // 0 - 

-> yp/cu) = 3.4 keV. 

(25) 


One should keep in mind that the results (24) and (25) are obtained by assum¬ 
ing the scalar vertex to be pointlike. The procedure which allows one to include 
the effects of finite-range scalar meson formfactors in a gauge-invariant way 
is well known [8,29,30]. As shown in detail in [28], these corrections are small 
in the case of the (f) —> yao/Zo decay. The reason for this is the following: the 
(f) as well as the a 0 //o are close to the KK threshold, and the loop integral is 
saturated by nonrelativistic values of the loop momentum, \k\ <C where 
rriK is the kaon mass. The range of the scalar formfactor is defined by the 
range of the force and, in the absence of pion exchange between kaons in the 
scalar sector, the latter is obviously larger than the kaon mass. On the con¬ 
trary, the mass of the p/ui is significantly smaller than 2 thk, and the typical 
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values of the momentum in the loop integral are not that small. Thus, for the 
decays S —> yp/co, one expects corrections due to the finite range of the scalar 
formfactor, which would reduce the pointlike result. 


4 Comment on the 77 decays of scalar mesons 


The transition a 0 //o —> 77, closely related to the class of the reactions S —> 
yV, probes the matrix element M(p 2 , q 2 ) in a kinematical regime quite differ¬ 
ent from the decays discussed above. 

The 77 decay can proceed via the quark loop mechanism. Nonrelativistic quark 
model estimates give [38] 

r„( 3 a) = ^r 77 ( 3 P 2 ) = 432a 2 (Q 2 ) 2 iih|7 (26) 

where R(r) is the radial part of the wave function, and M 0 is the mass of 
the P- wave quark-antiquark state, which, in the leading nonrelativistic ap¬ 
proximation, is supposed to be the same for all the members of the P- wave 
multiplet. 

The calculation of the squared charge factors ( Q 2 ) 2 yields for the isospin ratios 

T[(a 0 -*• 77 ) : (f 0 (nn) -> 77 ) : (f 0 (ss) -*• 77 )] = 9 : 25 : 2. (27) 

Therefore, one can try to estimate the decay width for the scalar / 0 (980) from 
the width of the tensor / 2 (1270), which is known to be a good nn state. The 
PDG [34] quotes 


T(/ 2 (1270) -»• 77) = 2.61 ± 0.30 keV, 


(28) 


that gives 3 

r(/„(980) ^ 77) = f (sjf)) r (A< 1270 ) - 77) = 4.5 keV, (29) 

Similar results were obtained in other computations of T(/ 0 (980) —> 77) based 
on the qq model of the scalar mesons [39,40]. 

In the qqqq picture, the predictions appear to be of the order of 0.3 keV for 
both the /o and the a 0 [41]. 

3 Although the final state contains two photons and, therefore, the matrix element 
scales as lo 2 , the phase space brings the factor of lu/M 2 , with M = 2 uj being the 
physical quarkonium mass in its rest frame. Therefore, the relation T oc a ; 3 holds. 
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The 77 decay can also proceed via the kaon loop mechanism (see, e.g., Ref. 

[42] ), with the matrix element given by Eq. (18) with rriy = 0. For a pointlike 
scalar with the mass of 980 MeV and g 2 s f 4-7T = 0.6 GeV 2 one obtains 

T(S -> 77 ) « 0.24 keV. (30) 

In line with the reasoning of the previous section, this value comes out as our 
prediction for the 77 decay of a molecule. However, in the kinematical regime 
of the S —> 77 transition, the momenta in the kaon loop are in the order of 
the kaon mass and, as in the case of the S —» 7 V transitions, one expects 
corrections due to the finite range of the form factor at the scalar vertex. 
Indeed an explicit calculation within a molecular model of the scalar mesons 

[43] yields results for the decay widths (0.20 keV for the / 0 and 0.78 keV for 
the do) that differ from our predictions, but are still in remarkable qualitative 
agreement with them given the simplicity of our approach. 

The result of Eq. (30) as well as the applied technique is very different from 
those in Refs. [39,44], where also the 77 width of scalar molecules was calcu¬ 
lated. The authors obtain T(f 0 (KK) —■> 77 ) = 0.6 keV [39] and 6 keV [44], In 
these references, similarly to the positronium 77 decay, the transition matrix 
element is taken proportional to the value of the KK wave function at the 
origin. Not only is this quantity model dependent (as reflected in the order of 
magnitude variation of the calculated widths), but also we suppose that the 
validity of such an approach is highly questionable for the considered decays: 
the range of the KK —» 77 transition operator is of the same order as that of 
the wave function. 

The experimental values for the 77 widths of scalars are [34] 

r 77 (/ 0 (980)) = 0.391°;}° keV, T 77 (a 0 (980)) = 0.30 ± 0.10 keV. (31) 

The estimates based on the nonrelativistic quark loop (29) are in clear dis¬ 
agreement with these data. Although relativistic corrections to the formula 
(26) evaluated in Ref. [47] reduce the ratio r 77 ( 3 P 0 )/r 77 ( 3 P 2 ) by a factor of 
2, this result is still much larger than the experimental values (31). On the 
other hand, the kaon loop mechanism estimate (30) is certainly compatible 
with them. Moreover, as shown recently in Ref. [45], the new data [46] for 
the reaction 77 —> 7 r + 7 r~ in the vicinity of the /o(980) resonance can be de¬ 
scribed with the kaon loop mechanism using the weight factor Ds(m 2 s ) which 
reproduces the .S-wave tht scattering data. 

Concluding, one can say that the existing data on the 77 widths of scalars 
seems to favor a molecular structure of the /o/oo mesons. However, one has 
to admit that, at present, no reliable estimation of the theoretical uncertainty 
involved in the value quoted in Eq. (30) can be given. 
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quark loops 3> meson loops 


mg MeV 


Fig. 5. Illustration of various kinematical regimes probed by the decays involving 
scalars. 

5 Discussion 


The results of the previous sections are summarised in Fig. 5, where, for the 
sake of transparency, the relative contributions of the quark loop and kaon 
loop mechanisms are displayed for various kinematical regimes probed by the 
radiative decays involving scalar mesons. 

The message disclosed by this figure is quite clear: the closer the mass of the 
vector meson is to the KK threshold, the larger is the contribution of the 
kaon loop mechanism to the decay amplitude. Now recall that the quark loop 
mechanism is of relevance only if the scalars indeed carry a significant quark 
component, whereas the kaon loop mechanism contributes in both cases, i.e., 
in the qq (or qqqq) as well as in the KK molecule scenario. Thus, in order to 
discriminate between these two scenarios, it is most promising to study those 
decays where the quark loop mechanism, if present, is significant. 

The data on cf radiative decays yield [34]: 

Br((f) —■> qao) — 7.6 • 10 -5 , T(0 —■> -ya 0 ) — 0.3 keV, 

Br(c/> -> 7 / 0 ) ^ 4.4 • 10" 4 , T(0 -> 7 / 0 ) ~ 1.9 keV, 

indicating that the kaon loop mechanism indeed dominates the radiative tran¬ 
sition ej) —> 7 S. This can be understood from the proximity of both the mass 
of the ej) and the mass of the scalar mesons to the KK threshold. The ques- 
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Table 1 
Radiative 

tion that remains to be addressed is, however, how much room is there for an 
additional quark component. We argue that we can get a handle on this quark 
component of the scalar structure when looking at the decays S —> yp/co, 
for there the kinematical situation is quite different. As seen from the esti¬ 
mates given by Eqs. (11) and (25), in those decays the quark loop and kaon 
loop mechanisms should yield contributions of the same order. Accordingly, 
for scalar mesons with a qq structure, where both mechanisms contribute, 
the radiative transition widths should be significantly larger than for KK 
molecules, where only the kaon loop mechanism can occur. Thus, the radia¬ 
tive decays ao/./o —> 7 p/u> appear to be a much more decisive testing ground 
for discriminating between models for the scalar mesons than the radiative <f 
decays. 

The estimates of the radiative decay widths of scalar mesons in various models 
are collected in Table 1. The numbers demonstrate that the qq component of 
the scalar mesons implies characteristic ratios for the radiative decays into 
isovector or isoscalar vector mesons, namely T(a 0 —> yp)/T(ao —> yco) ~ 
1/10 and T (/ 0 —> 7 p)/T (/ 0 —> yco) ~ 10. This is in strong contrast to the 
corresponding ratios for the meson loop contributions (driven by the kaon 
loops), where all transitions are predicted to be of the same order of magnitude 
so that those ratios should be in the order of 1 . 

The latter point means that the scalar radiative transition is a filtering re¬ 
action. The quark loop mechanism “senses” the qq flavour, with the ratios 
of rates for different isospin content given by Eq. (13). Thus, one is able to 
measure the qq content of a specific scalar meson produced in a specific re¬ 
action simply by measuring the ratio of decay rates (/ 0 —>► yp)/(fo —* y^>) or 
( a o -> 7^)/(«o 7 P)- 

Another advantage of the radiative scalar decays is related to the fact that 
the phase space available for the final state is not small, in contrast to the 
radiative decays of the <f meson. Thus, simultaneous studies of data on f 
radiative decays and scalar radiative decays could be useful in establishing 
such fundamental characteristics of the scalar mesons as their pole positions 


Decay mechanism 

Process 

Radiative width in keV 

quark loop 

ao IP 

14 


a 0 yu; 

125 


fo(nn) -> 7 p 

125 


/o(nn) -> yu 

14 

KK loop 

ao/fo ->• yp/u 

3 


S —> yptio transition in various models. 
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and coupling constants. Indeed, the kaon loop mechanism is dominant in the 
0 radiative decay. The transition matrix element in the kaon loop mechanism 
exhibits a rather peculiar dependence on the masses of the initial and final 
mesons. The photon emitted in the 0 radiative decay is relatively soft, u ~ 40 
MeV, so that the corresponding matrix element decreases very rapidly in the 
upper part of the scalar-mass range, from the KK threshold to the mass 
of 0 (see the first plot in Fig. 4). On the other hand, in the reaction S —> 
yoj/yp, the photon energy appears to be large, about 200 MeV, so that the 
matrix element exhibits a rather different pattern. In the quark loop model 
it is nearly constant. As far as the kaon loop mechanism is concerned, the 
corresponding matrix element decreases rapidly with the scalar invariant mass 
from the KK threshold downwards, but remains practically constant above 
the KK threshold (see the second plot in Fig. 4). In such a case it is then 
possible to analyse the upper part of the spectrum with much less uncertainty 
than in the 0 radiative decay. 

For the sake of completeness we would like to mention that the radiative decays 
of scalar mesons have been also studied within the Vector Dominance Model. 
Corresponding results can be found, e.g., in Ref. [48]. 


6 Summary 


We have demonstrated that for radiative decays of the scalar mesons /o(980) 
and Oo(980) to the vector mesons p and u, meson loops and quark loops lead 
to very different predictions for the ratio of the decays to p and cu, respectively. 
Specifically, it follows from our results that for objects with a significant com¬ 
ponent from a compact quark state both types of loops should be equally 
significant. On the other hand, for scalar mesons that are KK molecules only 
meson loops can contribute. The inferred estimates for the decay rates and, in 
particular, the ratios that follow for the two scenarios are so drastically differ¬ 
ent that it should be possible to discriminate between them once experimental 
information becomes available. 

We have also pointed out that the radiative decay rates involving the scalar 
mesons (such as 0 —> yS, S —> yV, and S —■> 77 ) exhibit a distinct hierarchy 
pattern for a compact as well as for a molecular structure of the scalars. This 
pattern can be likewise used to distinguish between the two scenario, and, as 
an ultimate goal, to define the admixture of the bare confined state in the 
wave function of the scalar mesons. It requires, however, that a detailed and 
consistent calculation of all those rates is performed within a particular model 
for the scalar mesons. In this context let us mention that the molecular picture 
of the scalar mesons has been already successfully tested for the decays 0 —> 7 S 
and S —» 77, for which experimental data are available. Note that there exist 
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calculations [49] which reproduce both the value of 2.61 keV (see Eq. (28)) 
for the 77 width of the / 2 ( 1270) and also the / 0 (980) data (see Eq. (31)), in 
contrast to results of the nonrelativistic quark model (29). Obviously, it would 
be important to perform calculations within the quark-model picture of the 
scalar mesons for the other decays discussed in this work. 

Thus, experimental data on the transitions o 0 //o —> 7 p/oj — especially when 
analysed together with the existing data on </> —■> 7 a 0 //o — will provide strong 
constraints on models for the structure of the scalar mesons and, therefore, 
are an important source of information towards a solution of the scalar-meson 
puzzle. 
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A Differential width for unstable particles 


For the calculation of observables, like the decay width, in addition to M, two 
more ingredients are relevant — namely, the propagator of the scalar meson 
Ds{p 2 ) and that of the vector meson Dy(q 2 ). The latter modifies the invariant 
mass spectrum of the final state and one can use the unitarity relation to 
introduce the spectral function py for the vector meson, 

(2tt) 3 f d$k(pv]Pu ■ ■ ■ t Pk) \D v (q 2 )W v \ 2 = --Im D v (q 2 ) =: py(q 2 ), (A.l) 

J 7T 

where the integral denotes the integration over the phase space of the decay 
products of the vector meson that emerged from the vertex ITV• Note that 
the spectral density is normalized as 

J pv(q 2 )dq 2 = l. (A.2) 

The c o meson is quite narrow, so that choosing a Breit-Wigner form for p u 
is appropriate. For the p meson either a Breit-Wigner form or the data from 
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e + e —> 7r + 7r directly can be used. 

The finite width of the scalars makes one to study the decay rates as a func¬ 
tion of the invariant mass of the decaying system. Consequently, in the total 
decay width, Ds{m 1 2 3 4 5 6 s ) appears as a weight factor, in the m 2 s integration. For 
this distribution of the scalar mesons one would either need to use parame- 
terizations given in the literature or refer to data from the same production 
reaction where the radiative decay is extracted from. 

Having this in mind we can straight forwardly generalize Eq. (2) to the case 
of unstable particles in the final and initial state: 

j Pv(q 2 ) • (A. 

As mentioned before, the transition matrix element M is the quantity of inter¬ 
est and we investigate it for various scenarios in the main text. The theoretical 
predictions for the corresponding two dimensional distributions can be easily 
generated for each szenario discussed in the main text. 

It will be quite demanding to observe the radiative decays of a 0 and f 0 experi¬ 
mentally. First of all one has to identify reactions that allow one to disentangle 
the isoscalar f 0 and the isovector a 0 . Possible reactions that isolate, e.g., the 
former state would be dd —> a+scalar and J/T —> ^/cu+scalar. Then the 
intermediate scalar states needs to be reconstructed from the four-vectors of 
the decay particles. 


d 2 Y 


dq 2 dm 2 s 


m s 

327T 


\M(m 2 s ,q 2 )D s (m 2 s )\ 2 1 - 


mi 
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